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We characterise the kinematic image of the constraint variety of a 2R dyad 
as a regular ruled quadric in a three-space that contains a “null quadrilateral”. 

Three prescribed poses determine, in general, two such quadrics. This allows 
us to modify a recent algorithm for the synthesis of 6R linkages in such a 
way that two consecutive revolute axes coincide, thus producing a 5R linkage. 

Using the classical geometry of twisted cubics on a quadric, we explain some 
of the peculiar properties of the the resulting synthesis procedure for 5R 
linkages. 
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1. Introduction 

Given three poses of a rigid body, there exists, in general, a unique closed kinematic loop 
of four revolute joints with one degree of freedom, a so-called Bennett or 4R linkage, 
one of whose links visits the three poses. Suh (1969) solved this synthesis problem 
by constructing two 2R dyads to the three given poses and proving that they can be 
combined to form a Bennett linkage. The three pose synthesis problem for 4R linkages 
was re-visited by Brunnthaler et al. (2005) in a dual quaternion setting. A generalisation 
of their approach led to a factorisation theory for “motion polynomials” (Hegediis et ah, 
2013a). It was used by Hegediis et al. (2014) to synthesise spatial 6R linkages - closed 
kinematic chains with six revolute joints - to four prescribed poses in general position. 
More precisely, the four poses determine two families of (not necessary real) rational 
coupler motions, each giving rise to nine families of overconstrained 6R linkages. It 
should also be mentioned that the construction is only capable of producing a special 
type of 6R linkages. 
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In this paper we intend to close the gap in this sequence of exact factorisation based 
synthesis procedures for closed loop linkages and present a construction of 5R linkages 
(Goldberg linkages) to “more than three but less than four poses”. Our synthesis al¬ 
gorithm is based on a characterisation of the image of open 2R chains under Study’s 
kinematic mapping and the observation that the coupler motion of a 5R linkage is also 
parameterised by a cubic motion polynomial (see Hegediis et ah, 2013b). This allows to 
treat 5R synthesis as a special case of above mentioned synthesis of 6 R linkages. 

We continue this text with a quick introduction to dual quaternions and kinematics 
in Section 2. In Section 3 we provide a geometric characterisation for the kinematic 
image of a 2R dyad. While necessary properties can be derived easily, the proof of their 
sufficiency requires more work and comprises the main part of this text. The synthesis 
of 5R linkage in Section 4 is then merely a corollary to previous results. Note however, 
that it is not obvious how to use the additional degrees of freedom in exact 5R synthesis. 

2. Preliminaries 

This article’s scene is the projectivised dual quaternion model of spatial kinematics. 
Here, we give a very brief introduction to this model for the purpose of settling our 
notation. More details will be introduced in the text as needed. For more thorough in¬ 
troductions to dual quaternions and there relations to kinematics we refer to (Klawitter, 
2015, Section 3) or (Selig, 2005, Section 11). Familiarity with Hegediis et al. (2013a) is 
recommended too. 

The dual quaternions, denoted by OH, form an associative algebra in where mul¬ 
tiplication of the base elements 1, i, j, k, e, ei, ek, and ek is defined by the rules 

= k^ = ijk = —1, ^ = 0, ie = d, je = ej, ke = ek. 

An element q G OBI may be written as q = p + ed with quaternions p,d G HI := (1, i, j, k) 
(angled brackets denote linear span). In this case the quaternions p and d are referred to 
as primal and dual part of q, respectively. The conjugate dual quaternion is q = p + sd 
and conjugation of quaternions is done by multiplying the coefficients of i, j, and k with 
— 1. It satisfies the rule qf = fq for any q,r G OBI. The dual quaternion norm is defined 

ll^ll = qq- We readily verify that it is a dual number, that is, an element of O := (1, e). 

We identify linearly dependent non-zero dual quaternions and thus arrive at the pro¬ 
jective space = P(M®). Writing [q] for the point in that is represented hy q G OBI, 
the Study quadrie is defined as 5 :={[<?] G P^: ||g|| G M}. With q = p + sd, the algebraic 
condition for [( 7 ] G S is pd + dp = t). 

Identifying P^ with the projective subspace generated by (1,ei,ej,ek), a point [q] = 
Ip -|- sd] G S with non-zero primal part acts on [x] G P^ via 

[x] !-)■ [{p -|- ed)x{p — ed)]. (1) 

The map (1) is the projective extension of a rigid body displacement in Composi¬ 
tion of displacements corresponds to dual quaternion multiplication. This isomorphism 
between SE(3), the group of rigid body displacements, and the projectivisation of the 
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group of dual quaternions of real norm and non-zero primal part provides a rich and 
solid algebraic and geometric basis for investigations in kinematics. 

3. A characterisation of 2R spaces 

We complement the Study quadric by the quadric AA = {[g]: ||g|| G eM} of points repre¬ 
sented by dual quaternions whose norm has vanishing primal part. It is a quadric of rank 
four. The set i? = {[( 7 ]: g G eH} of its singular points is a projective space of dimension 
three and it is contained in the Study quadric. Only the points of E are real points of 
M whence we are led to consider the complex extension P(C®) of . We call M the 
null cone and refer to the set E of its singular points as the exceptional three space. 

Given two non coplanar lines fi, £2 in Euclidean three space, we consider the set of all 
displacements obtained as composition of a rotation around £ 2 , followed by a rotation 
about £ 1 . Its kinematic image is known to lie in a three space (Selig, 2005, Section 11.4) 
which we call a 2R space. By construction, a 2R space contains the point [1] (the 
identity displacement) and intersects 5 in a doubly ruled regular quadric. The rulings are 
obtained by varying one revolute angle and fixing the other. In order to characterise 2R 
spaces among all three-dimensional subspaces of with these properties, we introduce 
the notion of null quadrilaterals. These are spatial quadrilaterals contained in both, S 
and Af. We call their edges null lines. 

Theorem 1. A three space P C P^ with [1] £ P is a 2R space if and only if it 

• intersects the Study quadric in a regular ruled quadric Q, 

• does not intersect the exceptional three space E, and 

• contains a null quadrilateral. 

The first and second item in Theorem 1 exclude exceptional cases with coplanar, 
complex, or “infinite” revolute axes, the latter corresponding to prismatic joints. The 
important point is existence of a null quadrilateral. We split the proof of Theorem 1 
into a series of lemmas. It will be finished by the end of this section. 

Lemma 1. The straight line [x] V [y] is contained in S if and only if xx = yy = 
xy + yx = 0 . 

We omit the straightforward computational proof but note that the left-hand sides of 
each of the three conditions in Lemma 1 are dual numbers. Hence, they give just six 
independent conditions on the real coefficients of x and y. 

Lemma 2. The conditions of Theorem 1 are necessary for 2R spaces. 

Proof. The constraint variety of a 2R chain can be parameterised as 

R{tip2) = {ti — hi){t2 — h2) (2) 
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with two dual quaternions hi, h 2 that satisfy hihi = h 2 h 2 = 1 , /ii + /ii = /i 2 + /12 = 0 
and hih 2 7 ^ The first condition ensures that the dual quaternions are suitably 

normalised and [/ii],[/i. 2 ] G S. The second condition means that hi and /12 describe 
half-turns (rotations through an angle of vr). The third condition ensures that the axes 
of these half turns are different. Equation (2) describes a composition of two rotations 
for all values of ti, ^2 in M U { 00 } with 00 corresponding to zero rotation angle. Even if 
their kinematic meaning is unclear, we also allow complex parameter values. Expanding 
( 2 ) yields 

R{tl,t2) = tlt2 — tih2 — t2hl + hih2. 

We see that the kinematic image of the 2R dyad lies in the three space spanned by [1], 
[hi], [h 2 ], [hih 2 ]- In a suitable projective reference with these points as base points, 
the surface parameterised by ( 2 ) is the quadric with equation xqxz — xiX 2 = 0 which is 
indeed regular and ruled. 

The intersection of P with the exceptional three space E is non empty if and only 
if the primal part of R{ti,t 2 ) vanishes for certain parameter values ti, t 2 - This can 
only happen if the primal parts of hi and /i 2 are linearly dependent over M but then 
the revolute axes are parallel and P, contrary to our assumption, is contained in S. 
The other possibility for P C S, intersecting revolute axes, has been excluded as well. 
Clearly, P is not contained in the null cone M either. We claim that the intersection of 
P and M consists of the four lines given by ti = ±i, t 2 = ±i. Indeed, they are null lines. 
Take for example x = (i — hi){t 2 — / 12 ) and, in view of Lemma 1, compute 

XX = (i - hi){t2 - h2){t2 - /i 2 )(i - hi) = {t2 - h2){t2 - /i 2 )(i - hi){i - hi) 

= {t2 — h2){t2 — h2){—y^ — i(/ii -l- hi) + hihi) = {t2 — h2){t2 — /i 2 )(~l — 0 -|- 1 ) = 0. 

This derivation uses the fact that {t 2 — h 2 ){t 2 — h 2 ) is a real number and thus commutes 
with all other factors. The cases ti = —i, ^2 = ii are similar so that we have verified all 
conditions of Theorem 1. □ 

The proof of sufficiency is more involved. We need two additional lemmas from pro¬ 
jective geometry which are formulated and proved in Appendix A. 

Lemma 3. A three spaee P that satisfies all eonditions of Theorem 1 is a 2R spaee. 

Proof. Denote the vertices of the null quadrilateral (in that order) by [rti], [xi], [ii 2 ]) and 
[X 2 ]. Intersecting the Study null quadrilateral with the tangent hyperplane of S at [1] 
yields points 

[mi] G [ui] V [xi], [ni] G [xi] V [U 2 ], [m 2 ] G [U 2 ] V [X 2 ], [ 712 ] G [X 2 ] V [xi] 

of a planar and non-degenerate quadrilateral (Eigure 1). The joins [mi] V [m 2 ] and 
[ni] V [ 71 , 2 ] are real lines and their real points correspond to rotations about fixed axes. 
Pick real rotation quaternions [hi] G [mi] V [m 2 ] and [h 2 ] G [ni] V [ 772 ]. The axes of 
these rotations are the only candidates for our 2R dyad. Hence, we have to show that 
either [hih 2 ] G P or [h 2 hi] G P. It is easy to see that this is the case if and only if 
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either [mini] G P or [nimi] G P. In fact, we will even show that one of these products 
equals ui. 

We claim that both Mi := [mi] V [mini] and M 2 '■= [mi] V [nimi] are null lines. In 
order to show this, we have to verify the conditions of Lemma 1: 

(mini)(mini) = mi(nini)lnr = 0, (nimi)(nimi) = ni(mimi)lli = 0, 

=0 =0 

mi(mini) + (mini)mi = mi(ni + ni)tni = (ni + ni)(mimi) = 0, 

eC =0 

mi(nimi) + (nimi)mi = (mimi)ni + ni(mimi) = 0. 

=0 =0 

Similarly, we see that also Ni := [ni] V [mini] and N 2 '■= [ni] V [nimi] are null lines. 
Thus, we are in the situation depicted in Figure 1 where we have three null quadrilaterals 
with respective vertices 

[ui], [t>i], [U 2 ],[V 2 ]; [mmi], [m 2 ni], [m 2 n 2 ], [min 2 ]; [mmi], [nim 2 ], [n 2 m 2 ], [n 2 mi]. 

The second and third quadrilateral are different because mi and ni do not commute 
(otherwise they would lie on the same line through [1] which contradicts the regularity 
of the quadric Q := P r\S). Our proof will be finished as soon as we have shown that 
the first and the second or the first and the third quadrilateral are equal. For this, it is 
sufficient to show that [ui] = [mini] or [ui] = [nimi]. 

At first, we argue that the primal part of [ui] equals the primal part of [mini] or of 
[nimi]. Because P does not intersect P, the projection on the primal part is a regular 
projectivity P —)■ [H] with centre E. We denote projected objects by a prime, that is, 
we write u[, u'l, m'l etc., for the primal parts of ui, ui, mi etc. The quadric Q is regular 
and ruled and so is its primal projection Q'. Hence, the point [m'^] G Q' is incident with 
precisely two lines contained in Q', say M[ and M^. But [m'l] V [u'l] is contained in Q'. 
Hence [u(] G M( or [u^] G M^. Similarly, [i;(] is also incident with one of the two lines 
contained in Q' and incident with [n'l]. Thus, [u'l] = [rn'in'i] or [i;(] = [n-irn'i]. 

Now we have to lift this result to the dual part and show that [ui] = [mini] or 
[ui] = [nimi]. The alternative being similar, we assume [u^] = [min'i]. This means 
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that we have two null quadrilaterals, one with vertices [tti], [ui], { 112 ]-, [^ 2 ] and one with 
vertices [mini], [m 2 ni], [m 2 n 2 ], [min 2 ] such that their primal projections are equal and 
corresponding sides intersect in the vertices [mi], [ni], [m 2 ], [n 2 ] of a planar quadrilateral. 
But then Lemma 5 in Appendix A asserts their equality and [ui] = [nimi] follows. □ 

4. Synthesis procedures 

In this section we discuss the exact synthesis of 5R linkages, based on our geometric 
characterisation of 2R spaces. The attentive reader will note that we neither provide a 
synthesis algorithm nor a concrete example. The reason for this (and one of the main 
points of this paper) is that the synthesis of 5R linkages can be treated as special case 
of a recently introduced synthesis procedure for 6R linkages (Hegediis et al., 2014). In 
this context, 2R spaces play a crucial role. 

At first, we re-prove a classical result of Suh (1969) on the synthesis of 2R dyads to 
three poses. Its proof is just a variant of our computations in the proof of Theorem 1. 

Theorem 2. Three poses [po] = [1]? \pi], [P 2 ] £ S that span a plane that is not tangent 
to S and does not interseet E are ineident with preeisely two 2R spaees. 

Proof. By Theorem 1, the sought 2R space P contains a null quadrilateral. We compute 
its vertices [ui], [ui], [U 2 ], [^ 2 ] as follows. The plane [po] V [pi] V [p 2 ] and the Study 
quadric S have a regular conic C in common. This conic intersects the null cone in two 
pairs mi, m 2 , and ni, n 2 of conjugate complex points that can be computed from the 
roots of a univariate quartic polynomial. Again, we denote projection on the primal 
part by a prime. We can determine the primal part of ui by solving the system 
Uiu'i = u'lm'i + m'lu'i = tt'in 2 + Uin 2 = 0. It consists of two linear and one quadratic 
equation. Hence, there are precisely two solutions in projective sense. Picking one of 
them, the primal parts of the remaining vertices can be computed unambiguously: Write, 
for example, v[ = u'^ + Xm[ and solve + n'lv'i = 0 for A. Finally, from [rt)^], [u'^], 
[u 2 \, [^ 2 ] the points [ui], [ui], [U 2 ], [^ 2 ] can be uniquely recovered as in the proof of 
Lemma 5. □ 

Remark 1. The two solutions in Theorem 2 both contain the conic C. This conic is 
the kinematic image of the coupler motion of a Bennett linkage (Hamann, 2011). Thus, 
as already mentioned by Suh (1969), the two 2R dyads of Theorem 2 can be combined 
to form a Bennett linkage. Moreover, the computation of Lemma 3 shows that the two 
2R spaces are [1] V [hi] V [/ 12 ] V [/ 11 / 12 ] and [1] V [hi] V [/i 2 ] V [h 2 hi] with suitable rotation 
quaternions hi, h 2 . In other words, they differ only by the ordering of their axes. This 
shows that the figures formed by opposite axes in a Bennett linkage are congruent. 

In a 5R linkage there are two possibilities to assign a coupler to a fixed base. The two 
choices can be distinguished by the degree of their relative motion with respect to the 
base. As has been shown recently by (Hegediis et al., 2013b, Theorem 6), the degrees 
of relative motions between non-adjacent links in the dual quaternion model of rigid 
body kinematics are two and three, respectively. In order to gain additional degrees of 
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freedoms over the synthesis of Bennett linkages, it is mandatory to use the link to the 
degree three motion as coupler. In this sense, we can say that the coupler motion of a 
5R linkage admits a rational cubic parameterisation (a cubic motion polynomial) in the 
dual quaternion model of SE(3). However, not every twisted cubic in the Study quadric 
gives rise to a 5R linkage. The necessary and sufficient condition is that the coupler 
motion can also be generated by a 2R dyad, that is, the twisted cubic lies in a 2R space. 

Using the factorisation theory for rational motions (Hegediis et ah, 2013a), a generic 
cubic motion polynomial C can be written, in six different ways, as C = {t — hi){t — 
h 2 ){t—hz) with linear motion polynomials t—hi, i G { 1 , 2 ,3}, that parameterise rotations 
about fixed axes. Suitable combinations of such factorisations give an overconstrained 
6 R linkage whose coupler follows the prescribed cubic motion. For details on how to 
compute cubic interpolants on quadrics, how to factor the resulting motions and how to 
pick factorisations suitable for linkage synthesis we refer to Hegediis et al. (2014). We 
just highlight the major differences to the case of 5R linkage synthesis. 

By a fundamental property of the underlying theory, the factorisations of a cubic 
motion polynomial are, at least in generic situations, in bijection with the permutations 
of the real quadratic factors of CC (Hegediis et ah, 2013a). Pairings for “admissible” 
kinematic chains correspond to permutations with different factors at the beginning or at 
the end. In the case of 5R synthesis, the situation is special. At least one factorisation 
is of the shape C = {t — hi){t — h 2 ){t — / 13 ) where either the axes of hi and /12 or 
of /i 2 and /i 3 are identical. Lets assume the first case. Then hi and /12 commute so 
that C = {t — h 2 ){t — hi){t — hs) is another factorisation, different in algebraic sense 
but identical in kinematic sense. It can be paired with just two of the remaining four 
factorisations to form a closed 5R linkage. If /i 2 and h^ have identical axis, above 
discussion can be repeated with C = {t — h^){t — h 2 ){t — hi) instead of C. 

For a more detailed discussion of some properties of 5R linkages in this context, we 
need some results of twisted cubics that can be found in classical texts (for example 
Salmon, 1882; Cayley, 1885). The intersection of a 2R space with the Study quadric 
is a regular ruled quadric Q. It carries two families of rulings. One family contains a 
ruling corresponding, via Study’s kinematic mapping, to all rotations about the axis at 
the base of the 2R dyad, the other family contains a ruling corresponding to the moving 
revolute axes. We may accordingly speak of a “first” and a “second” family of rulings. 
The twisted cubics on Q can be partitioned into two classes: Members of the first class 
intersect every ruling of the first family in two and every ruling of the second family 
in just one point. For cubics of second class, the situation is just the other way round. 
Finally, four points on Q determine two one-parametric families of twisted cubics, one 
of first class and one of second. A member of each family is uniquely determined by the 
choice of one further point on a ruling through one of the four given points. 

These considerations clearly show that four poses in general position cannot be reached 
by a 5R linkage unless their kinematic images span a 2R space. On the other hand, three 
poses in general position already determine, by Theorem 2, two 2R spaces but also a 
Bennett linkage. By above discussion, we have three degrees of freedom to find a cubic 
curve through three quadric points. In other words, the number of free degrees in 5R 
linkage synthesis exceed that of 4R linkages by three. Unfortunately, only restricted use 
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can be made of these. It is, for example, not possible to arbitrarily prescribe one further 
orientation: While there is a dual quaternion with prescribed primal part (orientation) in 
each 2R space, it will, in general, not lie on the Study quadric. An additional prescribed 
position is prevented for similar reasons. 

Finally, also the appearance of two types of factorisations with commuting factors on 
the left {hi, /i 2 ) or on the right (/i 2 , h^) can be explained by the two families of twisted 
cubics on Q. If the commuting factors are on the left, C'o(t) = {t — hi){t — / 12 ) is a 
quadratic parameterisation of the line [1] V [hi] and every point on that line corresponds 
to two parameter values. In particular, there exist a second parameter value to (besides 
00 ) such that [C'o(to)] = [!]• But then C{to) yields a second curve point (besides [1]) on 
[1] V [/is], that is, the cubic belongs to the second family. In similar manner we see that 
the cubic belongs to the first family, if /12 and h^ commute. 

5. Conclusion 

We characterised the kinematic image of 2R dyads and used it to specialise a recently 
developed synthesis procedure for 6R linkages in such a way that it yields 5R linkages. 
Even if sufficient degrees of freedom are available, it is not possible in general to find a 5R 
linkage that visits three and a half poses (three poses plus one position or one orientation) 
because of geometric obstructions. Thus, a natural selection of interpolation data for the 
envisaged coupler motion is not entirely clear. Hopefully, future research will provide us 
with more ideas in this direction. 

A. Auxiliary results 

Here, we prove two technical results. The formulation of Lemma 4 could be simplified 
but in its present form its application in the proof of Lemma 5 is apparent. 

Lemma 4. Given a three-space E C and four points [u'l], [u(], [u'f\, [v'f\ that span a 
three-space F C P'^ with E D F = 0, set Ui := [u'^] V E, Vi := [u(] V E, U 2 ■= V E, 
^2 := [^ 2 ] V E and consider four projections 

hl'-Ul^Vi, Vi\Vi^U2, tl2'-U2^V2, l'2:V2^Ui 

with respective centres [mi], [ni], [m 2 ], and [ 7 ^ 2 ]. If the centres span a plane L, the 
composition k := V 2 ° h 2 ° ° hi of the four projections is the identity. 

Proof. Note that dim Hi V Vi = 5 so that the projection /ii (and also vi, p, 2 , 1 ^ 2 ) from a 
suitable point is well defined. Take an arbitrary point [tti] G Ui but not in E and not in 
L and set [ui] := /7i([rti]), [U 2 ] ■= ni{[v 2 ]), [V 2 ] ■= h 2 {[u 2 ])- Because these points are not 
coplanar, we may take them, in above order, as base points of a projective coordinate 
system in G := [rti] V [ui] V [U 2 ] V [U 2 ] and complement them by four further base points in 
E and a suitable unit point to a projective coordinate system of P^. We select the unit 
point such that its projection into G from E gives ([mi] V [m 2 ]) n ([ni] V [ 712 ]) and denote 
induced projective coordinates in Ui, Vi, U 2 , V 2 , and G (obtained by dropping three 



or four zero coordinates) by writing the subspace as subscript. With this convention, 
the projection centres are [mi] = [1,1,0,0]g, [m] = [o,i,i,o]g, M = [o,o,i,i]g, 
[^2] = [1)0,0,1]g- Moreover, we denote the coordinate vector of an arbitrary point 
X G C/i by [xo,*];7i where the star abbreviates four fixed coordinates. We then have 

[xo,*]c/i ^ ^ [xo,*]u2 ^ [-xo,*]v2 ^ [xo,*]ui 


and the claim follows. □ 

Lemma 5. Let E, F C be non-interseeting three spaees, the latter spanned by 
points [tti], [u'l], [u'^, [U 2 ]. Consider further a regular quadrie Q eontaining E and four 
vertiees [mi], [ni], [m 2 ], [n 2 ] ^ Q of a planar quadrilateral. Then there exists a unique 
spatial quadrilateral with vertiees [rti], [ui], [U 2 ], [^ 2 ] that is eontained in Q, whose sides 
[tti] V [ui], [ui] V [U 2 ], [U 2 ] V [v 2 \, [^ 2 ] V [ui] are, in that order, ineident with [mi], [ni], 
[m 2 ], [n 2 ] and whose profeetions from E into F are [u'l], [v'i\, [rt^, and [v'f\. 

Proof. Denote a quadratic form associated to Q by oj. Given [u'^], [u^], [^ 2 ], [v'f\, we 
have to reconstruct [ui], [ui], [^ 2 ], [^ 2 ] subject to the constraints 

uj{ui,ui) = u(vi,Vi) = u:{ui,Vj) = t), i,je {1,2] (3) 

and 

[mi] G [«i] V [ui], [ni] G [ui] V [ii 2 ], [^^ 2 ] £ M V M, K] G [U 2 ] V [wi]. (4) 

If [iii] is given, we find [ui] by projecting [lii] from centre [mi] onto E V [ui]. Lemma 4 
tells us that we can find [U 2 ] and [^ 2 ] in similar manner such that (4) is satisfied. Then, 
some of the conditions in (3) become redundant and it is sufficient to consider only 

Oj{ui,Ul) = Uj(vi,Vl) = U){U 2 ,U 2 ) = 0j{y2,V2) = 0. (5) 

In a projective coordinate system with base points [rt'i], [u'l], [^ 2 ], [v'f\ G F and further 
base points in E the quadratic form oj (and the quadric Q) is described by a matrix of 
the shape 

U B 
B O 

where A, B and O are matrices of dimension 4x4, O is the zero matrix and B is regular. 
Now (5) gives rise to a linear system for the unknown coordinates of [lii], [ui], [U 2 ], and 
[^ 2 ]. The matrix of the linear system is equivalent to the matrix B by elementary row 
transformations, whence existence and uniqueness of a solution hinges solely on the 
regularity of B. □ 
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